Current rectification by asymmetric molecules: 
An ab initio study 
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We study current rectification effect in an asymmetric molecule HOOC-CeH4-(CH2) n sandwiched 
between two Aluminum electrodes using an ab initio nonequilibrium Green function method. The 
conductance of the system decreases exponentially with the increasing number n of CH2. The phe- 
nomenon of current rectification is observed such that a very small current appears at negative bias 
and a sharp negative differential resistance at a critical positive bias when n > 2. The rectification 
effect arises from the asymmetric structure of the molecule and the molecule-electrode couplings. A 
significant rectification ratio of ~38 can be achieved when n — 5. 
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PACS numbers: 85.65.+h, 73.63. -b, 36.40.-c 
I. INTRODUCTION 

Electronic transport properties of single molecule junc- 
tions have gained tremendous interest in recent years, 
since they may have a wide variety of important appli- 
cations in future electronic components such as tran- 
sistors, diodes and switche s 1 ' 2 ' 3 ' 4 ' 5 . With the advan- 
tages of experimental techniques, for example, scan- 
ning tunneling microscope and mechanically controllable 
break junctionals, measurement of the current through 
nanoscale systems is now allowed. Some interesting be- 
haviors, such as highly nonlinear I-V characteristics, neg- 
ative differential resistance (NDR) and electric switching 
behavior, are observed in various systems ^ 10 ' 11 . Mean- 
while, considerable amounts of theoretical work have 
been performed to study transport properties of molecu- 
lar devices. 

Rectification, one of the most important functions in 
a traditional electronic component, has also been sug- 
gested and observed in molecular devices. Some mecha- 
nisms for rectification phenomena in these devices have 
been suggeste d 12 ' 13 : 14 ! 15 . The first molecular rectifier 
was proposed by Aviram and Ratner in 1974 using D- 
a-A molecules, where D and A are, respectively, an elec- 
tron donor and an electron acceptor, and a is a covalent 
"sigma" bridge (insulator). In this model, the inelastic 
electron transfer is more favorable from A to D, rather 
than in the opposite direction, thus current rectification 
happens. However, electrical rectification was observed 
experimentally only recently in Langmuir-Blodgett (LB) 
multilayers or monolayers of the molecule C16H33Q- 
3CNQ sandwiched between metal electrode a 16 ' 17 : 18 . The 
underlying mechanism, which was first considered to be 
a possible implementation of the Aviram-Ratner mech- 
anism, is of somewhat different D-tt-A type. Although 
the molecule does show rectification, it behaves as an 
anisotropic insulator rather than a conductor because of 
the long alkane tail C16H33. In fact, the current is of 
order 10~ 17 A /molecule for the structure AI/C16H33Q- 
3CNQ/A1 and lCr 15 A/molecule for the structure Au/ 
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FIG. 1: (color online) Model structure of a two-probe system 
with a single molecule HOOC-C6H4-(CH2)n-S coupled to two 
Al(100) electrodes. 



Ci 6 H 33 Q-3CNQ/Au in Ref£. The current is too small 
for practical applications. To achieve effective rectifica- 
tion with a reasonably conductive molecule, one should 
avoid using molecules with long saturated groups that 
may lead the molecules to be prohibitively resistive. It 
seems attractive to use relatively short molecules with 
certain end groups to allow their selfassembly on a metal- 
lic electrode's surface. In Ref. 19, transport measure- 
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FIG. 2: The (I-V) curves for different number of CH 2 : (a) 
n = 0, (b) n = 1, (c) n = 2, (d) n = 3, (e) n = 5 and (f) 
n = 7. 



merit on the structure of two phenyl rings connected 
to electrodes with asymmetric alkane chains (molecule 
(CH 2 ) m -CioH 6 -(CH2)„) has been performed^. It sug- 
gests another new mechanism of molecular rectification, 
where a single electroactive unit is positioned asymmetri- 
cally with respect to electrodes and the highest occupied 
molecular orbital(HOMO) and lowest unoccupied molec- 
ular orbital(LUMO) are positioned asymmetrically with 
respect to the Fermi level. For majority of the applied 
voltage drops on the longer insulating barrier, the condi- 
tions for resonant tunneling through the active level are 
different for different directions of the applied voltages 
on the two opposite polarities. A maximal rectification 
ratio 35 is expected with longer alkane chains n = 10 and 
?7i = 2, the current magnitude is still not ideal, with or- 
der 10" 10 A Recently ideal rectification has been demon- 
strated in diblock oligomer diode molecules^ and single 
C59N molecules- . However, the current reporte d 14 ' 17 is 
still of the unsatisfactorily order 10 _10 A In this work, 
we perform a systematic study on the transport prop- 
erties of a single molecule HOOC-C6H4-(CH 2 ) n . Some 
of our findings arc interesting. First, a significant recti- 
fication ratio of ~38 is achieved at bias 1.2V for much 
shorter alkane tail (n — 5), and a feasible current of or- 
der 10~ 8 A is observed, two order larger than reported 
in Refs i 14 ' 20 . Second, the rectification effect persists in 
a much wider region than in previous investigation !] 14 ' 20 . 
The rectification does not degrade when n is varied. 

The paper is organized as follows. The the device 
model and computation method are presented in Sec. II; 
the results and discussions are given in Sec. Ill; and a 
conclusion is drawn in Sec. IV. 
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FIG. 3: The equilibrium conductance as a function of the 
number n. 



II. MODEL AND METHOD 

The model structure for our theoretical analysis is 
illustrated in Fig. 1. The single molecule HOOC- 
C6H4-(CH2)n with sulfur end is sandwiched between two 
Al(100) electrodes with finite cross section. The elec- 
trodes are chosen from the perfect Al crystal along the 
(100) direction, and the number of atoms in each atom 
layer is arranged in sequence 5,4,5,4 • • -. Four Al atomic 
layers (5, 4, 5, 4) are selected for the electrode cell. The 
single molecule HOOC-C6H4-(CH2)„ with sulfur end to- 
gether with four surface atom layers in the left elec- 
trode and three surface atom layers in the right elec- 
trode is chosen as the central scattering region, as indi- 
cated by two vertical white lines in Fig. 1. The struc- 
ture of every molecule HOOC-C6H4-(CH2)„ with sul- 
fur end is optimized via Hartree-Fock approximation in 
ArgusLab— i 26 ' 27 i 28 , and the energy convergence of 10" 10 
kcal I mol has been achieved. The molecule- lead distance 
is fixed to be a constant for all the systems. The distance 
of the H-Al(from the left terminal hydrogen atom to the 
surface of Al electrodes) is 1.5 A and the S-Al is 1.0 A. 

The calculations for transmission and Current- 
Voltage^- V) characteristics have been performed using a 
recently developed first-principles package TranSIESTA- 
C, which is based on the nonequilibrium Green's func- 
tion (NEGF) technique. The TranSIESTA-C, as imple- 
mented in the well tested SIESTA method, is capable of 
fully modeling self-consistently the electrical properties of 
nanoscale devices, which consist of an atomic scale sys- 
tem coupling with two semi-infinite electrodes as shown 
in Fig.l. Such a nanoscale device is divided into three 
parts: left and right electrodes, and a central scatter- 
ing region. In fact, the central region includes a portion 
of the semi-infinite electrodes. The external potential 
bias is included in the self-consistent calculation directly. 
Therefore, the effects of the bias voltage on the electronic 
structure of the system can be fully considered. De- 
tails of the method and relevant references can be found 
elsewher o 21 ' 22 ' 23 ' 24 . In our actual calculations, the con- 
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FIG. 4: The transmission T as a function of energy E for 
n = 0, 1 and 2. For comparison, the transmission for n = 2 
is magnified by 10. 

vergence criterion for the Hamiltonian, charge density, 
and band structure energy is 10 -4 and the atomic cores 
are described by normconserving persudopotentials. 



III. RESULTS AND DISCUSSIONS 

Firstly, the equilibrium conductance as a function of 
the number n of CH2 has been studied and is presented 
in Fig. 2. We find that the conductance decreases ex- 
ponentially as the number n of CH2 increases. The con- 
ductance is 0.15G , 0.043G , 0.0043G , 0.00026G if the 
number of CH2 n = 0,1,2,3 (Go = 2e 2 /h is the con- 
ductance quanta). The similar exponential decrease of 
conductance with the length has been found for carbon 
chains^. 

The I-V curves are plotted in Fig. 3 (the cases n = 4, 6 
are not shown in this work). As shown in Fig. 3, the 
current takes a large value of order /j,A both at the neg- 
ative bias and at the positive bias, When n — and 1. 
When n > 2, interesting phenomenon are observed. The 
current through the molecule structure at the negative 
bias is increasing slowly irrespective of the bias voltage 
increase, as compared to the quick increase of the cur- 
rent with the increasing bias voltage in the positive bias 
case. When the bias is positive, the current reaches its 
maximum value at a bias, which is different for different 
n systems. These are typical current rectification and 
NDR phenomena. A significant rectification ratio ~38 is 
found for the molecule system with the number of CH2 
n = 5 , and the current in the negative bias side is neg- 
ligibly small. We note that, for every configuration, a 
strong NDR can be observed as the bias is positive. The 
NDR in our system may find potential applications as 
switches. 

What are the underlying mechanisms for the current 
rectification and NDR in our molecule system? Why the 
current rectification can be observed only when the num- 
ber of CH2 exceeds 1? To answer these two questions, 
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FIG. 5: (color online) The LUMO and HOMO distribution 
on the molecule for n — 0,1 and 5 respectively: (a) LUMO for 
n = 0, (b)HOMO for n = 0,(c) LUMO for n = 1, (d)HOMO 
for n = l,(e) LUMO for n = 5 and (f)HOMO for n = 5. 
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FIG. 6: Transmission spectra, T(E,Vb), for the n = 5 case. 
The two crossing dotted lines indicate the bias window. 



one need to know detailed information of the transmis- 
sion spectra of such a system. For current rectification 
to occur, two essential conditions should be satisfied that 
the transmission resonance must be reasonably sharp and 
exists only at one vicinity side of the Fermi energy^. 
The transmission spectra T(E) are plotted in Fig. 4 for 
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n = 0, 1 and 2 cases. We can see from Fig. 4 that the 
resonance is rather wide for n = case and has a long 
tail near the Fermi energy. A direct comparison reveals 
that the resonance in the case n = 1 is wider than that 
in the case n = 0, but moves toward right. However, 
the resonance becomes much more narrow and can be 
seen only at the right side of the Fermi energy for the 
case n — 2. The reason is that, when n = 0, the sul- 
fur atom is directly attached to the benzene, thus the 
orbitals of the lone electron pair of the sulfur are over- 
lapped with that of the ir electrons of the ring, resulting 
in molecular orbitals distributed almost over the ring and 
the sulfur. Such an observation is confirmed in Fig. 5(a), 
where both the LUMO and HOMO extend almost over 
the whole molecule. Since the sulfur is directly coupled 
to the electrode, the molecular orbitals are broadened. 
In the n — 1 case, the sulfur is separated from the ring 
by a single insulating -CH2- group. However, the separa- 
tion is not large enough to prevent elongated p orbitals of 
sulfur and sp3 hybridization of carbon from direct over- 
lapping, as seen in Fig. 5(6) where the LUMO and HOMO 
still distribute almost over the whole molecule. As a re- 
sult, the electronic level is still much broadened. Things 
changes when n — 2. Two insulating groups move the 
sulfur away from the ring by 4.09A, so that the direct 
overlap between the sulfur and the ring wave functions is 
weakened, and the LUMO and HOMO do not distribute 
on the left (CH2)2 (which is not shown). That's why 
rectification can be observed only for n > 2. 

To probe clearly the origin of rectification and NDR, 
we take the molecule HOOC-C6H4-(CH2)s as an exam- 
ple. When the molecule is coupled two Al electrodes, the 
calculation shows that the LUMO and Homo energies are 
0.394eU, — 2.896eU, respectively, consistent with the in- 
dications in the transmission spectra in equilibrium in 
Fig. 6. We can see from Fig. 5(c) that the LUMO dis- 
tributes over the HOOC and the carbon atoms in ben- 
zene ring, and the HOMO mainly on the sulfur atom. 
The electron has an accessible LUMO near the Fermi 
energy only of one electrode, the left electrode in the 
present work. This asymmetry leads to a significant cur- 
rent when the applied bias is positive and a negligibly 
small current in the negative case. 



The calculated transmission spectra T(E, Vb) are pre- 
sented in Fig. 6 for different applied bias voltages from 
— 1.5V to +1.5U with an interval 0.1V and a doubled 
interval within [-0.6,0.6]V. The dashed line represents 
the transmission spectrum in equilibrium. The cur- 
rent is obtained using the Landauer-Biittiker formula: 
1 = i*T T (^' V b) dE - Tne bias window contributing to 
the current integral is indicated by two dotted lines in 
Fig. 6. The transmission resonance in equilibrium ap- 
pears on the right side of the Fermi energy, and its width 
is narrow, about 0.13 eV. With the increase of the pos- 
itive bias, the resonance moves toward the bias window 
and makes contributions to the current. As a result, the 
current increases rapidly and reaches its maximum value 
at the bias 1.2 V. However, when the applied bias is 
negative, the resonance moves far away from the Fermi 
energy as the bias increases, and finally is pushed out of 
the bias window. It results in a negligibly small current. 
This is a typical phenomenon of current rectification. A 
notable rectification ratio of 38 can be observed at the 
bias 1.2 V. When the positive bias continues to increase, 
the transmission within the bias window is degraded in 
magnitude and width, leading to an NDR phenomenon. 



IV. CONCLUSION 

A system of HOOC-C6H4-(CH2) n is studied systemati- 
cally using a first principles method based on NEGF. The 
equilibrium conductance decreases exponentially with 
the increase of n. Strong rectification effect is observed 
when n > 2 and a significant rectification of ~38 is 
achieved from n = 5 for bias <1.2V. Irrespective of the 
number of CH2, a strong NDR can be seen under positive 
bias. A detailed analysis on the origin of rectification is 
given from the distribution of the HOMO and LUMO 
and from the transmission spectra at various bias. 
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